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Abstract 

We consider topological quantum memories for a general class of abelian anyon models defined on spin 
lattices. These are non-universal for quantum computation when restricting to topological operations alone, 
such as braiding and fusion. The effects of additional non-topological operations, such as spin measurements, 
are studied. These are shown to allow universal quantum computation, while still utilizing topological 
protection. Our work gives an insight into the relation between abelian models and their non-abclian 
counterparts. 
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1 Introduction 

Anyons are quasiparticles that can exist on two-dimensional systems [1,2]. Their 
non-local behaviour is an interesting topic for foundational research in quantum 
physics [3,4]. Importantly, one may employ anyons for fault-tolerant quantum com- 
putation [5,6,7]. The so-called non-abelian anyon models receive most attention, 
since they possess straightforward means to store, protect and manipulate quan- 
tum information. Even so, proposals have been made to realize universal quantum 
computation using an abelian model, the toric code [8,9,10]. We extend the latter 
scheme to a wider class of abelian models, proving universality and fault-tolerance. 
This provides fresh insight into the power of abelian anyons in relation to their 
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non-abelian counterparts. Physical realizations of abelian anyons are simpler than 
those for non-abelian anyons, giving an experimental motivation for our work. 

We employ Kitaev's quantum double models [5], expressed as stabilizer codes 
[11]. By analysing the outcome of stabilizer measurements, known as the syndrome, 
the type of error may be determined and corrected. The probability of the syndrome 
giving misleading results is suppressed by certain parameters of the code. We specif- 
ically consider models based upon the cyclic group of d elements, Z^, realized on a 
square lattice with a d-level spin on each edge. The elements g G are used to 
label basis states of d-dimensional spins, for which the generalised Pauli operators 
are defined, 

(1) E \9 + ^){9\, a^=Y.^'\9){9\- 

Here lo = e*^'^/'^. These satisfy the commutation relation a^a^ = uja^a^. The 
eigenbasis and eigenvalues of the operation are obvious. The eigenstates of 
are those of the Fourier transform basis, 

(2) \~g) = ^Y.^"'\h)^ 

with corresponding eigenvalues lo~^ . To rotate between these two bases, we use the 
Fourier transform, 

(3) F= \~9){9\ = ^=j E ^"'\h){9\, 

g&Za ^"g,/ieZd 

which has the properties F"^ \ g) = \ —g), = and -F^ = /. 

The stabilizers of the quantum double models are defined on the four spins 
around each vertex, v, and plaquette, p, 

(4) Aiv) = al^a^alal, B{p) = a^alaW 

where the numbering proceeds clockwise from the top- most spin (Fig. 1). These 
have eigenvalues lo^ = e^'^'^sl'^ for each g G Z^. No anyon is associated with a 
vertex or plaquette if A{v) \ ijj) = B{p) \ = | ■0). An anyon Cg is associated with 
a vertex if A{v) | V') = I V^)) ^-iid an anyon nig is associated with a plaquette if 
B{p) = u;^ I ■0)- The anyonic vacuum corresponds to the stabilizer space of the 
code, and forms the degenerate ground state space of the model's Hamiltonian, 
which may be expressed as, 

(5) H = -Y,A[v)-Y,B[p)+h.z.. 

V p 

This assigns energy to anyons, making them localised quasiparticle excitations. 

The creation and movement of the anyons may be achieved by and op- 
erations. The operation on spins 1 or 2 of a vertex, or a^^ on 3 or 4, creates 
an Cg charge at that vertex and an c-g on the other vertex shared by the spin. 
Similarly, a on spins 2 or 3 of a plaquette, or a a^^ on 1 or 4, creates an nig flux 
on that plaquette and an ni^g on the other plaquette shared by the spin. Particles 
can be moved and braided using corresponding strings and loops of the and 
operations. The commutation relations of these give a phase ui^^ when an eg anyon 
is moved clockwise around an nih. The phase is obtained for an anticlock- 
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wise braiding. Superpositions of anyon states may be moved using entangling gates 
[12,13] or potential wells [9], allowing more complex braiding behaviour. 




Fig. 1. The models are realized on a square lattice with a d-level spin on each edge. The numbering of 
spins around a plaquette, p, and vertex, v, is shown. Anyons Sg reside on vertices. These may be created 
in pairs by single spin operations, and moved away from each other by chains of these operations along a 
string of spins. Anyons mg reside on plaquettes, and may be created and moved similarly. 



2 Encoding of Qudits 

Logical qudits are stored on the lattice using states of anyon occupancies. The 
qudits come in two types, denoted v and p, stored on pairs of vertices and plaquettes, 
respectively. We consider a basis for these qudits labelled | j)^/p for j = . . . d — 1. 
These are defined in terms of the anyon occupancies of two vertices, vi and V2, or 
two plaquettes, pi and p2, as follows, 

(6) I j)v = I I , I j)p = I rnj)p^ \ m.j)^^ . 

The stabilizers are not enforced on the vertices and plaquettes in which the qudits 
are stored. This opens so-called holes in the code [10], bringing the logical states into 
the stabilizer space. The Hamiltonian (5) must be modified accordingly, with the 
corresponding plaquette and vertex terms removed. When the anyon quasiparticles 
are moved, e.g. for braiding, the holes must move with them. 

Pauli operators for the qudits are defined in a similar way to the lattice spins, 

(7) X = ^|j + l)(i|, Z = Y,^^\j){j\. 

j 3 

Certain logical operations follow naturally from the properties of the anyon model. 
Measurements in the logical Z basis may be performed by measuring the occupancies 
of the vertices and plaquettes. Logical Z and X rotations may be performed by 
braiding and fusing with ancillary anyon pairs, respectively. The Xg operation may 
be performed on the v qudits by creating an ei and e_i, placing the former on v\ 
and the latter on V2- Similarly, we may perform an Xp operation by placing fluxes 
on plaquettes. The Zg operation may be performed on the v qudits by creating 
an mi and braiding the mi around si and then annihilating the flux pair. 

Similarly, we may perform Zp by braiding a charge pair around pi. Braiding the 
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contents of the vi of one qudit around those of the pi of another will implement a 
controlled-Z gate or its inverse, depending upon the chirality of the operation. 

Note from (6) that the logical states | O)^^^ will take the form of Bell pairs. This 
anyonic entanglement is a natural consequence of the conservation law requiring that 
anyons are created in pairs. It has been shown that anyonic entanglement is truly 
non-local, since it may be used to violate Bell's inequalities [15]. In our case, the 
entanglement allows the qudits to be stored non-locally, providing fault-tolerance. 

3 Non-topological Operations and Universality 

The logical operations implemented by braiding and fusion do not form a powerful 
gate set. Further logical operations may be achieved by measuring the spins of 
the lattice. This allows logical qudits stored on neighbouring plaquettes or vertices 
to be measured in the X basis. It also allows the preparation of a wide range 
of logical states. We demonstrate how these may be used as ancillae in protocols 
implementing various single qudit gates. When d is prime, we prove the universality 
of the resulting gate set. 

Consider two logical qudits, one v-type and one p-type, stored on the vertices 
and plaquettes neighbouring a lattice spin, i. The vertex vi is taken to be that for 
which i is labelled 1 or 2 and the plaquette pi is that for which i is 2 or 3 (Fig. 2). 
The logical X operations then take the following simple form, 

(8) Xy = af , Xp = erf. 

Any projector onto states of i may be decomposed as a sum of the spin Pauli 
operators, and therefore the logical qudit Pauli operators. Single spin measurements 
on i therefore correspond to measurements of the logical qubits. In general, these 
will be in an entangled basis. The exceptions are measurements of af and af, which 
provide X basis measurements of the v- and p-type qudits, respectively. 
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Fig. 2. Measurements of a single spin i effect the surrounding plaquettes and vertices. A ii-type qudit 
may be stored in the two vertices and a p-type qudit in the two plaquettes. The single spin measurements 
correspond to either single qudit measurements, or entangling two qudit measurements. 

The X measurements may be used to implement a Fourier transform, F, on a,v 
type qudit in arbitrary state | ijj)^ = Ylj I J%- This requires an additional p-type 
qudit in state | O)^ = | k)^ l^fd-, which can be prepared probabilistically by Xp 
measurements. Entangling these together with a controUed-Z yields the state, 

(9) A(Z) I V^) J 0)^ = 5] u^^^c, I j). I = 4^ E I 0. ^^v^ I ^^v')- 

j,k ^ I 
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Measuring the v qudit in the X basis then teleports the state to the p-type qudit 
while implementing the Fourier transform F. This is up to a Pauli correction, 
according to the measurement outcome. A corresponding process can be used if the 
initial state is on a p-type qudit. The circuit for the process is given in Fig. 3(a). 
The Fourier transform allows the application of a controUed-X between two f-type 
qudits, or two p-type qudits. This is done by the application of F'^ to the target 
qudit, followed by the controlled- Z, and finally followed by F. 

Single spin measurements may be used to prepare logical ancilla states as follows. 
First a lattice spin, i, is measured in a basis including the state — ^ ^^'^^ I j)- 
The corresponding projector may be expressed, 

(10) I </>)(</> I = 1 ^ S'f>.+.-<t>,) I j + J,), (j I = ^2 E e'^^^+^-'f'^^u:-'J'xlxl 

j,k j,k,l 

This results in a superposition of the neighbouring v and p qudit states, whose co- 
efficients depend on those of Applying xj, to the p qudit and projecting the 
its state onto | 0)^, will leave the v qudit in the state 



Vd 

A Fourier transform using the above method will then transform this into the p 
qudit 

3 

Setting (/)o = and = — Ylik=i j-i for j = ^■■■d — 1, this becomes 

(11) \0), = Y.^''^\^)v^ 

3 

which is the required ancilla state. 

The ancilla states may be used to perform single qubit rotations of the form, 

(12) f^.,p(^) = E ^^'^ I 0- 1 , 

i 

to p-type qudits in arbitrary state | ■0)^ = Ylt \ k)^. To do this, the inverse of a 
controlled-X is applied between | ip)^ and an ancilla in state | 6)^, with the latter as 
the target. The resulting state is 

(13) E ^'''^^ I I i - = E ^^''^'^^ I I j% 

3,k 

Measuring the second qudit in the Z basis and obtaining the outcome j = means 
the phase gate has been implemented correctly. Otherwise the process may be 
repeated with the same qudits until the right result is obtained, changing the Oj^s 
to correct the erroneous rotations. Each attempt succeeds with a probability of 0.5, 
so success will come within a small number of steps. The circuit for the process is 
given in Fig. 3(a). 

These phase gates, along with the Fourier transform, are able to perform a large 
number of single qudit rotations. For the case that d is prime it is proven that 
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arbitrary single qudit rotations are possible [14]. With the entangling gates, this 
provides universal quantum computation. 
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Fig. 3. The circuits to implement (a) a Fourier transform, and (b) the rotation U{9). Both of these are 
implemented on a i;-type qudit in arbitrary state | i/i). Corresponding circuits exist for p-type qudits. 



4 Fault-tolerance 

We consider two forms of error in the model: perturbations in the Hamiltonian 
and thermal error. The former only affect the encoded information if they lift 
the degeneracy. Due to the non-local encoding of the qudits, such perturbations 
must consist of many body operations that either loop around the holes or stretch 
between them. These become increasingly unlikely as the distance between the holes 
is increased. Any other perturbations will not commute with the Hamiltonian, and 
so their effects are efficiently suppressed by the gap. 

Thermal errors applied to the spin lattice correspond to the propagation of stray 
anyons. The creation of these costs energy, and sothe gap natually protects against 
them. Once created, the movement of stray anyons does not cost further energy. The 
fault-tolerance of the encoding then comes from the underlying stabilizer code [6], 
since measuring the syndrome allows stray anyons to be detected and annihilated. 
If the anyons propagated far enough the annihilation will not correct the error, but 
consolidate it. For example, consider a f-type qudit stored in two holes, vi and 
V2, d vertices apart. If the syndrome is measured and a single e_i anyon is found, 
its antiparticle must reside in one of the holes. It will be assumed that it resides 
in the hole closest to the e_i, since this is the most probable case, and so the e_i 
will be moved there to annihilate the ei and correct the error. However, if the ei 
actually resided in the hole furthest from the e_i, this action would not correct the 
error, but actually form a logical X or X'^ operation on the qudit. Such a process 
requires errors to occur on d/2 spins, and so has probability 0{j)'^l'^\ where p is the 
probability of a single spin error. Moving the holes far apart therefore exponentially 
suppresses X errors. A corresponding argument holds for p-type qudits. 

The occurrence of Z errors requires the braiding of flux anyons around holes. 
Since the boundary of a hole consists of four spins, only two spin errors are required 
to fool the stabilizer and cause logical errors. This is because such errors form 
half of a loop around a hole, and attempts at error correction may complete this 
loop by mistake. The solution to this is to use a repetition code, which emerges 
naturally from the anyon models. To do this for 't;-type qudits, two additional holes 
are introduced. The first is on a vertex neighboring wi, denoted v'^- The shared spin 
between these is projected onto the state | 0), corresponding to the identity element 
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of Z^. To do this, the spin may be measured in the basis. If the outcome is \g)^ 
the operation A^^{y']) may be appHed, where the sign depends upon the orientation 
of the Hnk. This then maps the spin state to | 0) without creating flux anyons. The 
second additional hole is on v^-, neighbouring f2, and the shared spin is similarly 
projected to | 0). The projections cause the logical qudit to be stored in all four 
holes, and so it becomes further delocalised. Repeated applications of the procedure 
allow the qudit to be stored in 2N holes, arranged in two rows. The application of a 
logical Z requires flux anyons to braid around entire rows, rather than single holes. 
For this to occur in error, an anyon pair must propagate around at least half of the 
loop, which has probability 0{p'^/'^). The addition of holes therefore exponentially 
suppresses Z errors. The probablility of X errors, however, increases since there 
are more holes that stray anyons may fall into. This increase is polynomial, and so 
may easily be compensated for by increasing the distance between the rows. 

This scheme is called a repetition code since the X basis measurements can be 
performed using any two vertices from different rows, and with loops around single 
holes affected only the X-basis measurements using that hole. Measurement of the 
syndrome is then equivalent to majority voting with these. A corresponding scheme 
can be used to store p-type qudits on more holes, though the projections in that 
case are to the state | O) using B^3{p[) operations. 

The preparation of the logical ancilla states requires the holes to be neighbouring, 
and so prone to errors. The preparation should therefore occur on quickly, on 
timescales for which the probability of error is small. In addition to this, errors in 
preparation can be suppressed by using purification schemes [16]. 

5 Application to Other Models 

The methods developed above have a straightforward generalisation to all other 
quantum double models. These are similarly realized on spin lattices, with stabiliz- 
ers defined on each plaquette and vertex [5] . Flux quasiparticles are associated with 
non-trivial eigenvalues of the plaquette stabilizers and charge quasiparticles are as- 
sociated with those of the vertices. For the abelian quantum double models, those 
based upon an abelian group, these quasiparticles are gauge invariant and may be 
directly identified with anyons. For the non- abelian models the quasiparticles are 
gauge dependent the identification is not so straightforward. Even so, it is known 
how the gauge invariant anyons may be formed from these [12]. 

Logical information may be stored by not enforcing stabilizers, and so creat- 
ing holes. For the abelian models, the logical information corresponds the anyon 
states stored within the holes. For the non-abelian models, the encoding could be 
expressed either in terms of the anyons, or the underlying flux and charge quasipar- 
ticles [17]. The separation between the holes provides fault-tolerance, with errors 
suppressed by distance. The usage of additional holes will also provide protection for 
the abelian models and in some cases for the non-abelian, but the full generalization 
is not known. 

For any quantum double models, projections onto lattice spins may be decom- 
posed into creation operators for quasiparticle pairs. Spin measurements may then 
be used to measure logical states stored on neighboring plaquettes and vertices, and 
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to create ancilla states. These will then increase the computational power of the 
model from that possible by braiding and fusion alone. Note that, for non-abclian 
models, it may not be possible to measure v-type and p-type qudits separately, 
corresponding to the X measurements above. It is likely that only entangling mea- 
surements between a v- and a p-type qudit will be possible in general. 

The application of the methods beyond quantum double models is speculative. 
We must therefore use general principles to guide us. In any physical realization 
of anyons, there are two types of entanglement to consider: that of the underlying 
physical medium and that of the anyons. Any controlled reduction in the former 
seems to yield increases in the latter. For example, the state with vacuum every- 
where has no anyonic entanglement, but requires a maximally entangled state of 
the physical medium. Conversely, a separable state of the medium results in indef- 
inite anyon occupancies across the system. The corresponding superpositions lead 
to highly entangled states of anyons. Localised measurements on the medium will 
lower its entanglement. Anyonic entanglement should then arise as a consequence, 
which could be used to prepare logical ancilla states and increase computational 
power. 

The method can also be applied to non-anyonic stabilizer codes. Consider the 
following code, defined on six spin- 1/2 particles. 

(14) 5i = afal, S2 = ^3 = a|af , ^4 = afaf, 

A qubit may be stored in holes by not enforcing stabilizers Si and 5*2. The logical 
state I 0) is identified with a -|-1 eigenvalue for both Si and >S'2. The state | 1) is 
associated with a —1 eigenvalue. The logical Pauli operators are Z = Sioi S2 and 
X = (jf. Alternatively, a qubit may be stored by not enforcing ^2 and ^3. The 
Pauli operators in this case are Z = S2OTS3 and X = cr^crf. In the former case 
X is provided by a single spin operation, and so the information stored in this 
way is easily accessible. This allows the possibility of creating states by single spin 
measurements, as with anyons. In the latter case the X, like the Z's, is a two spin 
operation. This protects the stored information, since single spin operations are not 
sufficient to cause a logical error. 

With anyonic codes, the separation of holes means that logical qudits can be 
given arbitrary protection. In the above example the protection is limited to single 
spin errors. This limitation is also found when information is stored in holes of other 
stabilizer codes, such as the Steane and Shor codes. Of course this is not suprising, 
since these codes were designed to have information stored in the stabilizer space, 
rather than in holes. It would be interesting to know whether topological codes are 
the only ones for which holes provide arbitrary protection, or whether others are 
possible. 

6 Experimental Realizations 

The use of simple abelian models means that our work is in tune with current ex- 
perimental proposals, as well as actual experiments. One promising possibility is 
Josephson-j unctions, for which the means to produce the Hamiltonians for D{Zd) 
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models has been well studied [18,19], along with other models such as the non- 
abelian D^S^). Polar molecules have also been considered as a means to produce 
the Hamiltonians of quantum double models, with a complete toolbox of interactions 
proposed [20]. Topological states may also be produced without a Hamiltonian. Ex- 
periments with these have been performed, demonstrating anyonic braiding [21,22]. 
These systems, though small, could also be used to demonstrate the preparation of 
ancilla state described here. 

7 Conclusion 

We have studied a class of anyon models, realizable on spin lattice. These are 
abelian, and so not expected to possess any serious computational power. We show 
that the addition of single spin measurements increases the computational power 
significantly. This is because they may be expressed in terms of the logical states 
as measurements in entangled bases. This additional resource provably provides a 
universal gate set for many of these models, and likely for most others as well. We 
also comment on the fault-tolerance, naturally provided by the Hamiltonian and 
stabilizers of the topological models on which the schemes are based. This work 
also gives an insight into the relationship between abelian and non- abelian models. 
The repetition code used, for example, bears striking similarities to the encoding 
found in the non-abelian charges of the D^S^). Elaborations on these points may 
be found in further works by the authors [23,24]. 
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